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ABSTRACT 

The second term has been obtained in the asymptotic series 
for the second (longitudinal) adiabatic invariant of charged 
particle motion in a static magnetic field. This correction to the 
lowest order invariant has two sources: the correction to the 

lowest order magnetic moment and the integrated effect of the 
guiding-center drift across the field lines. The second term is 
found to vanish at the mirror points; therefore during its motion 
between mirror reflections, the guiding center deviates from the 
surface on which the lowest order invariant is constant and 
intersects this surface at reflection. 
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I . INTRODUCTION 

This paper contains a systematic derivation of the next 
term beyond the lowest order for the second adiabatic invariant 
of a charged particle in a magnetic field. A charged particle in an 
electromagnetic field possesses up to three approximate invariants 
of its motion. The first is the Alfven invariant 1 or magnetic moment 
^^/2moB, whereigj^is the component of particle momentum nw perpen- 
dicular to the magnetic field B.. If there is an electric field, p^is 
the momentum in the reference frame in which the electric field 
vanishes. Magnetic moment invariance requires that fields vary 
slowly compared to the particle gyration period and gradually 
compared to the gyration radius. The second invariant 2 ’ 3 is an 
invariant of the guiding-center equations of motion, which are 
equations obtained by averaging over the particle gyration about 
the magnetic field line. The second invariant is therefore, also 
an invariant of the particle motion, from which the guiding-center 
equations are derived. This second invariant exists when the 
guiding-center motion along a field line is nearly periodic; the 
invariant is j^ds , where ^ is the component of guiding-center 
momentum parallel to B>, and the integral extends over a period of 
the motion in s, which is distance along the line. The third 
invariant 3 is an invariant of equations of motion resulting when 
the guiding-center equations are averaged over the periodic motion 
along the field line (and therefore is an invariant of the particle 
motion also); it exists when these doubly averaged equations have 
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nearly periodic solutions, and this occurs when the surfaces of 
constant second invariant are topologically cylinders. The third 
invariant $ is the magnetic flux threading a second invariant 
cylinder. Its invariance is trivially true in static fields and 
becomes of significance only in time -dependent fields. 

The three invariants described above are really only the 
lowest orders of three asymptotic series of the form 
constant s Mo + eM 1 + e 2 M 2 4- , 

constant a» Jo + eJi + e 2 J2 + , (l) 

constant =■ $0 4 - + e §2 + , 

where Mo is ^/ 2 moB, Jo is <j>^pds, and $0 is the magnetic flux threading 
through a surface of constant Jo. The expansion parameter e is the 
mass-to-charge ratio of the particle. The invariance of Mo, Jo, and 
$0 can be surmised (and demonstrated) by rather physical methods, 
while to obtain higher order terms in each series may require deep 
insight or a systematic method. In the present paper we use a 
systematic method due to Kruskal^ to obtain Ji for static magnetic 
fields. Although the second invariant exists whenever the motion 
along a field line is nearly periodic, we confine ourselves to the 
case where the motion is oscillatory between two mirrors. The case 
where the particle traverses a closed field line always in the same 
sense must be treated separately in the systematic derivation. 

However, a direct derivation of Ji to be given at the end of the 
paper shows that the result is the same as for the oscillatory case. 

The proof of the invariance of Jo in reference 3 is valid for the 
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closed field-line case as well as for the oscillatory case. 

II. THE SYSTEMATIC METHOD 

The systematic method presupposes coupled equations of motion 
of the form: 

dx/dT = f (x, e) , (2) 

where x. has a finite number N of components and £ possesses a power 
series expansion in the small parameter e. Furthermore, it is 
required that all solutions of the system dx/dT = £(x., 0) traverse 
closed trajectories in £ space. The equation of motion of a charged 
particle can be put into this form with e equal to m/e. The method 
shows how to obtain a transformation from the N variables to another 
set (:z, 0) which have the property that the equations of motion are 

dz/dT = eh^, e), 

d0/dT = uj(z_, e) , ( 3 ) 

the important point being that 0 does not appear on the right sides, 
and that x s x(z, 0, e) is periodic in 0. The vector z has N-l 

'Vtn AM 

components. The transformation and the new functions h^ and <jj are 
obtained as series in e. 

If in addition the equations (2) are canonical, with x^ the 
vector (j 3 , q) , then there exists the adiabatic invariant 

constant = 0> e ) * ^ ^ 

where p and q are the canonical momentum and position, and 



where the integral is over a period of 0. The invariant is 
obtained as a series in e and, although a function of the z 

/Wv 

variables, may be rewritten in terms of the original x variables 
by inverting the transformation. To have an invariant, it is not 
really necessary that (2) themselves be canonical-only that they 
be transformable into canonical equations. 

The transformation x to (jz, $) is not made directly, but 
for convenience by way of intermediate variables (y, u) , where y 

/W 

is any vector constant on the closed, lowest -order (s =r: 0) 
solutions of (2) , and u is an angle-like variable specifying the 
position around such closed curves. 

The theory as described so far will produce only one adiabatic 
invariant series (the magnetic moment), but to any order in e 
desired. The second (and third) adiabatic invariants are obtained 
from "reduced" equations of motion as follows: it can be shown 

that the Poisson bracket [0, M] equals 1, which means that 0 and M 
can be used as conjugate variables in a canonical transformation 
from (p, q) to (p 1 , M; q’, 0) where p 1 and q 1 have each one less 

''Wv -w\ 'VWV 'Wx vvs 

component than p and q. (M is the sum of the magnetic moment serie 

^ am 

to the order in e to which one is working). The new Hamiltonian 
^'(p 1 , q 1 , M, e) is independent of 0 since K equals — dfc-'/d?), ancl 

*** AW 

M is zero, being the first invariant. Thus the reduced system of 


equations is: 
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•« _ d*U l (.R l > flj, M 3 e ) 
q i a dp! 

( 5 ) 


Pj - - aJ., M, e ) 

dq! 


If these again have the property that all solutions are closed 
in (p 1 , q') space when e is zero, then a second invariant exists 
in terms of new z variables, which we will call z 1 : 


J (z.' ) - fp' (z‘ , 0' ) 


(■&.' , 0 ' ) 
90 ' 


d0‘ 


( 6 ) 


New intermediate variables , u ! ) will generally also be used. 
Finally a second reduction can be performed to obtain the third 
invariant. The variables and transformations are illustrated in 
Fig. 1. Each transformation is carried out as a series in the 
expansion parameter e. 

Although (p 1 , q 1 , M, 0) are shown by the arrow in Fig. 1 
as coming from (p, q) s they really come from the entire first 

/v*. 

line, since (y, u) and (z, 0) must be found in order to determine 

A/%* ^ 

M. A similar statement holds from the 2nd invariant level to 


the 3rd. 
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With so many transf ormations and variables, tricks for 
reducing the labor are very welcome. There are several such 
shortcuts. For one thing, 0 is never needed, since it is simpler 
to use u as the integration variable in (4)- Likewise in (6), u‘ 
is to be used. Secondly, if the dz/dT equations of motion (3) 
have all solutions periodic when e (or another parameter) is zero, 
as for Eq. (2), or if they can be transformed so that this is so, 
the theory following Eq 0 (2) may be repeated. For the charged 
particle, the equations of motion (3) for are themselves of the 
prescribed form (2) with the needed periodicity when e = 0, after 
a trivial rescaling of the independent variable. Thus (^' , u 1 ) 
and (zj , 0') can be obtained from the variables rather than by 
the pathway of (p 1 , q f )« Furthermore, it is proved in reference (4) 

/yvv AM 

that : 



0 ') 


. La 1 > .-$' ) 

d0' 


d 0 1 = j>p(z' 


0 ', 0 ) 


• da. ( &' ?„#' ? 0 ) d 0> 

S0' 


(7) 


that is, the original canonical variables may be used. Thus the 

9L’ 

pathway of transformations actually followed in the present calculation 
is illustrated in Fig > 2. The transformation is between z and 
(y 1 , u 1 ) and not between (z, 0) and (y 1 , un- 


reduced canonical variables (jd 1 , 


) , and need not be found. The 




(y' , u')<— >(z' , 0') --- 

|AW ^ 


level of 1st 
invariant 

level of 2nd 
invariant 


Figure 2: The transformations and variables actually used in this calculation. 
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Here (jr , vj) are the position and velocity of the particle. The 
x variables can be regarded as either (p, q) or (r, v) , since 
the particle equations of motion can be written in the form of (2) 
with the required properties using either set. 

In (7) the right side appears to be a function of 0 but really 
is not, since the left side is not; the 0 dependence of the right 
side actually must disappear. This fact gives another shortcut -- 
namely, that any particular value of 0 desired can be used in the 
functions p(z', 0 1 , 0) and q(z ! , 0', 0). The choice 0 =* 0 greatly 
reduces the algebra. 

A fourth way to simplify the calculation is to choose the y 
variables cleverly. The y vector is required only to be a constant 

•w* 

of the lowest-order motion, and the (zj 0) follow uniquely from the 
(y, u) . There clearly are infinitely many suitable choices for they 

/wv rVVV 

variables, since any function of a given y is also a constant. The 
guiding principle is to choose for the components of ^ quantities 
which are both simple and constant to as high an order in e as 
possible, even though the theory only requires y to be a constant of 
the lowest-order motion. For example, one component of y 1 is much 
better chosen as Mo + e Mi rather than just Mo alone. Similarly, 
one component of y is much better chosen as the approximate guiding- 
center position than as the particle position, since the guiding 
center does not move as rapidly as the particle. 

Although this systematic method of calculating invariant series 
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in principle requires only routine labor, in practice the amount 
of algebra gets out of hand rather rapidly unless the above short- 
cuts are used, along with physical intuition in the choice of y. 
After this long but necessary discussion of the method, we now 
begin the actual calculation. 


III. THE CHARGED PARTICLE 
The particle equation of motion can be written as 


dr/dt a ev, 

dv/dT =* v X B(r) , (8) 

where x is (r, v) and t is t/e. When e is zero, r is constant and 

IVW -V-s Aw (VV\ 

the second equation (8) is the equation of motion in a uniform 

magnetic field B^. The solution has a constant projection along 

and is (harmonically) periodic in its two components perpendicular 

to J3. Thus the trajectory is periodic in space. Let L, M, and 

be three orthogonal unit vectors, with L^ parallel to M and need 
not be specified further. The y variables will be denoted by (p, a, T]) 

<vv* *** 

and defined as 


Z a z + e X * i$(£) /B h) > 

a = |vj_| , 

“H - k(jE) : 2* 

The choice of r itself instead of the guiding center was tried for 
p but led to a more difficult calculation. The y components T| and a 
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are the components of particle velocity parallel and perpendicular 
to the field line at the particle position. The u variable is 
chosen as 


u = 4: arc tan ^ 

2n v • N(r) 

Av. Vv 'A* * 


do) 


The inverse transformation can be obtained to any order desired from 


v - %(*;) + a [M(rJ sin u + N(r) cos u] , 


(ID 


r ^ p+-ecr 


j^)[N(r) sin u -M(r) cos u] 


where S3 *» 2rru. The equations of motion for ^ (du/dT is not needed) 
are : 


dCT/dT *-(7]e/a) vv:tfL(r) , 
dT)/ dT » evv:vL(r), 

A ^ A A A ^ A ^ 

dp/dT = eT)L(r) + ef_ il\ z L X (L -0L) + a 2 [ (M sin u 4 FJ cos u) X 

~ B l ~ ~ ~ ** 

/ A ^ A ^ A V A ^ A ^ A ^ A r> , 

( (M sin u + N cos u) * ) 4 (N sin u — M cos u) (M sin u + N cos u) • 

\ VA ' ** J V Ma \v* ' 5 j 

A ^ A ^ A A ^ y ^ ^ ^ yA\, 

4 T]a[M sin u + N cos u) X (L • VL) + L Xy(M sin u 4 N cos u) *VLJ ] V (12) 

1 aw ' V* ^ /v>^ y/w /w y 

where the : notation means contraction first of the two inner vectors 

A A A 

and then of the two outer ones. For example in (13b) MN : VL means 

i • a -v) Lj. 

The components of will be denoted by the Greek capitals of the 
corresponding y components - i.e., by (P., E, H) . By use of the 

'Vvk 

recursion and periodicity relations (B19 — B28) in reference (4) , we 



obtain 


P = p + 0(e 2 ) , 


( 13a) 


V — rT 4> - { ~ A ^ XV A /S <\ A A ^ 

- jrf - \i “T ( ~ M cos u + N sin u + M) L:VL + T|of — 4 MM sin u cos u + 

B *** *"* /Vu A* ' “ W AA 

A A A A A A a"\ 

^(MN + NM) sin 2 .u + ^ NN sin u cos t(] : VJL j- + 0(e 2 ) , ( 13b) 

r A A A A A A A 

H = T]+ e i — Tb( — M cos u + N sin "o + M) L : VL - a 2 !" — 'k MM sin u cos 

~7 s L ,Wv *»* ^ >* Afc 

B(j>) 


+i(MN + NM) sin 2 u + i NN sin u cos ul :X7Lf+0(e 2 ), ( 13 c) 

c VM' c *w*». J Arf* J 

0 = Jo ■'• (13d) 


The expression for 0 is not needed, since u will be used as the 

integration variable; (13d) is included only to show that 0-0 

when u = 0, a fact that will be used later. The € 2 terms of 

(j?, £, H) have also been calculated, but are too lengthy to include 

here. Note that choosing to be the guiding center has made the e 

contribution to J? vanish, and this will afford much simplification 

later. 

The equations of motion for the ^variables (d0/dT is not needed) 

are : 


-Jp = HL(P) + -I-, H 2 L X (L * VL) + §|z L X VB + L (? MN “ 1^) : VL 
G dT ^ B (P ) ^ V 2B ^ 2B ^ ^ 


+ B H 2L(ML^L) + 0(e 2 ) , 
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1 d£ HE 


6H- 


s dT ' 2 V ~ + B 


[B7 . I M . } £L ) _ ^(ml : 7 £)?*£] 


+ £ SlS (fcV* [L(MN + NM) : 7L] + BV • £(& : V7l) 

DO''. ^ 


B 2 


-B7 • 4£_L3zi)} t | H|i { g :V £ )v .£ + 0 


U 2 ), 


(14b) 


1 dH Zrj 
e dT = 2 V 


L - 


„ times the order e terms of — ^ + 0(e 2 ) 
H € dT ' ' 


(14c) 


Since st is time t, the left sides are the guiding -center 
velocity dP/dt, etc. The component of dP/dt perpendicular to the 
field is the sum of the usual gradient-B and line -curvature drifts, 
while the parallel motion has both zero-order and first-order 
components. The equations for d£/dt and dH/dt give 


d S 2 + H 2 
dt 2 


0(e 2 ) 


(15) 


which is energy conservation, while those for d£/dt and dEJ/dt give 


d 

dt 



(16) 


which is conservation of the lowest -order magnetic moment Mo. 

The so-called !, guiding-center equations* 1 , which are usual in 
numerical integrations designed to follow the guiding center, are 
the set (see for example reference 5) 
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S - H i(£) * if?) i * * w k 


^(l; + n ) = constant, 


( 17 ) 



constant . 


This is a hybrid set, in that only some of the terms of order € are 
retained in dP/dt and none of them in dZ/dt. It is simpler than 
the complete set (14) and its use has validity, as will be 
discussed at the end of the paper. 

At this point the next term of the magnetic moment series can 
be calculated. We will not give details, but just the result. 


Mo 4- eM, = — -Sjjfr [(NM L + £(NM + KN) :VL] 

2BIPJ dB #***» “Ki a»» * U< u,' am - 1 


«* IT ^ A A /N 

%-ML:7L 


( 18 ) 


Since the highest order term we will need in the magnetic moment 
series is , M will stand henceforth for Mo 4- *Mi , so that dM(zJ/dt 
is of order e 2 . When expression (13) is written in terms of (r,^v) 
it agrees with Eq. (28) of reference (6). 

We now proceed to the reduced system and the second invariant. 

If we set e - 0 in (14) we have 


dP/dt = HL(P) 
ds/dt - - ^HEV'L , 
aH/dt -= 


( 19 ) 
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Figure 3 


z 



The projection of the lowest-order trajectory into the £ 


- H 

H plane. 
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The first equation says that the guiding -center moves along a 
field line without deviating from it, the drifts having vanished 
for e a 0. The second and third equations show how the parallel 
and perpendicular energies interchange in the usual mirror manner. 
The motion is the five-dimensional space (P, E, H) is a closed 
loop; the projection of the loop into the E — H plane is double, 
as illustrated in Fig. 3» 

The^ 1 variables are to be chosen as constants of the motion 
along this path. One component of y 1 is naturally taken as the 
energy K = (E 2 4- H 2 )/2. Another component is the magnetic moment, 
and in accord with the technique of choosing the y variables as 
constant as possible, we use M * E 2 /2B(P) 4- eM-j , where Mi is 
given in (l8). The final two components of are a(P) and b(P) , 
where a and (3 are two functions of position, constant on each 
field line, such that the vector potential is and B is 

(See references 7 and 8). Since the lowest-order motion is strictly 
on a field line, a and (3 are sutiable for y' components. The angle 
variable o' we define as the fraction of its total longitudinal 
oscillation period T the particle has completed: 


U'(S, H, P) 

.Vn 


T(E,H,P) 


Po(E,H,P)* 2 ~t 2 






L(P') -dP 1 

T 


2B(P) 


eM^S.H.P) 


1 


b(p' 

•VA 


( 20 ) 
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where the path of integration is along the field line on which J? 
is located, and Pq is a zero of the denominator. For the oscil- 
latory case we are considering the denominator will vanish at 
two points Po and P 1 . For definiteness, choose Po as the one at 

Ah ah v* 

which B is directed towards rather than away from the other, as 

Vi 

in Fig. 4* The period T is twice the integral from JPq to F ^ . The 
positive sign in the denominator is to be used when integrating 
from ^Po to J?! and the negative sign on the return. A point on a 
field line has two values of u' whose sum in unity. It should be 
noted that the denominator in the integral (20) is not exactly 
the parallel guiding-center velocity at j? 1 ; it is to lowest order 
only. From (14) the parallel guiding -center velocity is 


£ £,-jj£ • H + (|mn - tm) :VL + SjjS (ML:VL) . 


This can be solved for H, which can then be substituted into 


2[K - (Mo + eM , ) B] . The result is 

2[K - (Mo + «M,)B] -jv^ 2 + e Vjl ^gi (NM -MN):(7 l] + 0(e 2 ). (21) 

Therefore the denominator of (20) is not exactly v^. When v^ is not 
too small, the square root of (21) may be expanded to give 

2^(K — MB)^= v n +2eMo(NM -MN):tfL + 0(e 2 ). 

|{ A*** 
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In the special case where J3>7* B - 0, the difference 

1 * 

between the parallel guiding-center velocity and ±2^(K —MB) 5 ' 

^ A ^ ^ A ^ 

does vanish. The reason is that (NM — MN) :VL equals (L/B) 

The vanishing of B*VxB is the necessary and sufficient condition 
for the existence of a family of surfaces orthogonal to the B field. 

■W\ 

Simplification sometimes appears in adiabatic theory when this 
condition is met. For two other cases, see reference 7, pages 30 
and 70 . 

In Fig. 4, the guiding center is shown reversing its parallel 
velocity at P- (i.e. at u 1 - |) because by (21) the denominator of 
(20) vanishes up to terms of order e 2 when v - 0^ and P 1 is defined 
as a zero of the denominator. The order e 2 difference is invisible 
to the order to which we are working. 

The equations of motion (14) in terms of (y 1 , u 1 ) are 

or = J? • Vat . 

3 = p -V3, 

A"* 

M = 0( e 2 ) , (22) 

K = 0(e 2 ) 



where the dots mean d/dt and where is given by (14a); H, and J? 

are to be expressed in terms of (a : y g , M, K, o'), a procedure that is 
in principle possible given the form of the field. 
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By using the systematic method for finding the z' variables, 
we obtain 


zi ' 


a +T^M) J‘ o du[J’ o d U "P(u") *17a(o") 


- P(u>?a(u)] + 0(e 2 ) , 


z 2 ' = p + T^Kp) J^dtCj’du" P(u")-7S(u") -P(u)-^(u)] + 0( e 2 ) . 
z 4 ' = K + 0(e 2 ) , 


(23) 


z 3 ' = M 4 - 0(e 2 ) = Mo + eM, + 0(e 2 ). 


1 

The integral P d’tf'P(u" )•?(*( u n ) is (to the order needed) the average of 

J 0 ^ /V - 

ot over the unperturbed path and will be denoted <<i>, and similarly 
for 3. The equations of motion for the 2 ^ variables have not been 
obtained from (14) because they are not needed in calculating the 
second invariant. 

The second invariant (since p ^ mv -b eA(r)) is 


iLai) 






d £(a' 1 , u' .i. p ) , 

&u' 


(24) 


where (7) has been used and the integration variable switched to u 
from 0 1 . The integral (24) is independent of 0. The range of 
integration for either 0' or u' is 0 to 1. The right side of (24) 



must be evaluated through order e so as to get the term of the 

series Jo 4- eJ-j 4- . The contribution to order 1/e in (24) 

vanishes because it is 1/e times the magnetic flux through the 
zero-order trajectory. But the zero-order trajectory is along a 
field line and back along the same path, thus enclosing no flux. 
The contribution to order 1 in (24) will be Jo, so we need up to 
and including order e. 

To get , u 1 , 0) correct through order e 2 and v(z', u’, 0) 

correct through order e is straightforward and tedious. Great 
simplification occurs when we set 0=0. When 0=0, u = 0 also, 
as shown by (l3d). From (l3a-c) we then have that P =5 p , S ^ a, 
and H = T- In fact according to the general theory these are 
exact relations (to all orders in e) when 0=0. Then from (ll) 


p eai(x ) p mil) , 0( ^ 

= 1 ~ B(r) = 5 B(P) + P H ° U } 



HL(r) + SN(r) = HL(P) + £N(P) 


B(P) B(P) ' 


The next step is to express (JP , 2, H) in terms of (z f , o'). In 
principle it is possible to invert the definitions of the (y 1 , o') 

fVVV 

to get (P, £, H) in terms of (y 1 , u')- Let 


P = R(a, p, M, K, u') 

be the formula for P obtained by inversion. Then P(z', o') is 

/Vw 

obtained by substitution of (23) into the funtion R and then Taylor 


expanding : 



( 27 ) 
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Hz', u') = R(zi ' , Z 2 ' , z 3 ' , z 4 ' , u' ) - T(<d>u' - f‘ U d) 

/*>'***'>** (J Q 

M ill ’ » Z2 ! , Zxl a Z4 1 . U') - T (<g>U ' — [° P) 

« “ 0 * ^ ' ’ 


where the T, d, and 3 are to be expressed as functions of (zj , o'). 

Consider the contribution of the vector potential to the integral 
in ( 24 ): 


r du* 

J n 


A (p ( z 


_ e ^ (P) + -®Lriji.jLrj 

~ V ~ ; ' E 2 B J du'«- 


e£M 


Vsf (M-VM - MM.-VB ) , 

+^7- - - B 


(28) 


where P^ 1 , \j') is to be replaced from (27). Since the integral 

(28) is needed through order e 2 , it seems that the term of P(z' , u') 

would be needed. As the following analysis will show, neither it 

A A A A 

nor (e 2 £ 2 /B 2 ) (M*VM — MM*VB/B) contributes and all of the order e 2 
contribution to the integral arises from the products of first order 
terms. The path of integration is at constant z* . Varying u 1 at 
constant z‘ does not make r(z', o') follow a field line unless 

yVl " 'W- 

o[r(z' , u')] and g[ < r(z', u')] are independent of u' ; but they do depend 
on u' - In fact 


Qf[r(z' , u' )] 


zi 


— T(<q , >u i 



iMi), 

B 


(29) 


and 



Z2 1 -T(<g>u' 



eHl-73(R) • 

B 


(30) 
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The ring of constant z 1 might look as in Fig. 4* The 
vector-potential integral in (28) is the magnetic flux passing 
through this ring, and that is the negative of the double integral 
JJdadg over the ring, because A. is aV0 . In an a - 0 plane the ring 
might look as in Fig. 5- The double integral is the area of this 
ring, which by a little geometry is 

JJdadg = y o du’{[e(z', u') -g(z’ , 0)]|jj, -[a(z', u') - a(z' , 0)]^j,}, 

(31) 


where a and ^ are given in (29) and (30) as functions of (z' , u 1 )- 
The difference cf(z ! > U 1 ) — a(z' , 0) is order e 5 since the zero order 
of each is zy 1 . Moreover 9a(z ! , u')/Bu' is order e. Similar 
statements hold for g(z', u‘)- These facts are also clear from Figs. 

4 and 5> where the deviation a(z', o') — a?(z^ , 0) of the ring from 
the U* = 0 field line is due to the drifts, which are of order e. To 
summarize, (31) vanishes through order e while its e 2 part comes only 
from the products of e terms. 

/S A 

We will next show that only the ej( M*Va)/B and eE(M *\7(3 ) /B terms 
of (29) and (30) contribute to (30; products of e terms coning from 
the a and g parts all cancel. Substitution of (29) and (30) into (3l) 
gives 

Jj’dadg =2 J 0 dU ' C (<3 >u ' ~ f U g)(<a > ~ a) — (<ar>u’ - a)(^P> -p)] 


+ eT j* du'f (<6>u' - f u 'a)-2- 
o J o L -'o B 
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+ f J 

p 1 du ' 
0 

j3<a> _ a) - 

(<8> 

-s) 

eT 

2 c 

fdu' 
• 0 

r«&> -a) (adZa. ' 

)u'«*0 

- (<e> - e)l 

The last integral vanishes because J 

' («*> _ ot) : 

0 


( 32 ) 


The second integral in (32) can be integrated by parts, and after 
some cancellation one finds 


r i 


r 

Jo 


du' ■ -J*J*dadp = eTj^du'l [ (<g> -p)M-tfa - («*> -a)M*V3] 


+ (<0xau'> -<«xpu'>) 


i2 1 


+ 2 J du , C<a>J' u du ,, p(u") - <6>J U du"of(u") ] 


+ j 1 du'[g(u' )J U du"a(u") - a(u')j' U du ,, g(u")]' 


(33) 


We proceed to show that the terms with T 2 give zero. In the second 
T 2 integral interchange the order of the u' and u" integration to 
get 

1 1 • «*.• 

J o du'<Qt>J^du"p(u") = <QOj o du"p(u") (l - u") = <«xg> - <«><8 u> (34) 


and similarly for the other part of the integral. Thus the T 2 terms 
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in (33) become 


T 2 {<gXau'> - <aX$u'> - ^ du'J U 'du"[a(u')e(u") - p(u' )a(u")]J. 


(35) 


Further progress depends on a theorem: Given two functions f(u) and 

g(u) such that f(l-u) equals f(u) and g(l-u) equals g ( 1 ) - g(u) . 

1 1 

Then f dug(u)f(u) = ^g(l)J* duf(u). In (35) the first term is 


1 . 1 


if 0 ^J ® v ' * Clearly, o' is a g function. Also, a i s an f function; 

this is because 1-u 1 is the same point on the zero-order field line 

• • 

as o', and the drift velocity^ which produces a a nd 8^ is independent 

of the sign of the parallel guiding-center velocity. Application 

1 . . * • 

of the theorem now gives i Qfu'scau^ = 4<a>> and similarly <gu'> = ^<8>- 

0 

And so the first two terms of ( 35 ) cancel. In the last term we have 

P du ’g (u‘ ) P U du M a(u 1 ') , where ft (u f ) is an f function, and P U du n <*(u M ) 

J 0 J 0 0 0 

is a g function of its upper limit. Application of the theorem yields 


J du '8 (u 1 ) J^du M a(u") = |<a><8 > and similarly J du ' ot( u 1 ) J U du"|3 (u n ) = 

^<8><a>. Thus all the T 2 terms disappear and the e contribution of 
the vector potential to j(z')/m is 


eJ ACS(S'- U ')3 Su > 

0 

(<a> -a)M(R)*yp(R)]. 


U ' )du« = tJ du'^y [ (<6> -8)M(R)*V«(R) - 


(36) 


In this expression jt means R^( zi ' , z% 1 , Z3 * , Z4 1 , u* ) and S is 
[ 2z3 1 B(Rj]^ to the order needed. The a and g are to be expressed 


as functions of z 1 and u 1 . 
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We now turn to the evaluation of the integral in (24) containing 
the velocity. Substitute r and v from (25) into (24) • 



where (P , H, £) are to be expressed in terms of (z 1 , o'). Equation 

/Vv 

(27) gives P(z' 5 u')> H(z' , o') and S(z' , o’) can be obtained by 
inverting the definitions of the y 1 variables and then expressing 
y 1 in terms of (z 1 , u 1 )* 

The lowest order of (37) is found to be the usual longitudinal 
invariant 


— Ls^J= pb 2 i [ z4 . 


b(r )] 1 


where the plus sign is to be used for 0£u and the minus sign for 
' 4 1 ) because the parallel velocity, which to love st order is H, 
changes sign at the turning point where o' is The plu s -or -minus 

sign is essential; without it Jo would vanish because L • SRV has 
opposite signs but the same magnitude at u' and 1-u 1 . Because of it, 
many terms in Ji will vanish due to the symmetry properties of their 
integrands . 

1 

The e part of (37) is not difficult to calculate. Many terms 
vanish due to the symmetry of the integrand. A sample term is 
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A * A. /v 

2 Z 3 1 M:VN which vanishes because M*7N is the same at u 1 

0 OU *** / '* v ^ 

and l~u‘, but dR/du ! has opposite signs. The result we find for 

the contribution of the velocity to J/m is 

/VvCz 1 , „■) U ' W = J-’du'(±)2*C**' -z,'B(R)]iL(E)- 


-Tj>'|(<*>V C±2^(Z4 ' 

+ (<B>u‘ -;“‘e) -z,'B )K- ^ 

+ (<d> -d)N(R)'-^- (2 z 3 'B)^ 

/W A*A QZ 1 

+ (<B> -B)N(R)-g^V (2 z 3 'B)^ 

- e ;’d 0 '( ± ) 2 i(z*' -z,'B)*(^)4[L^-:TO +L • 


(39) 


The last integral vanishes because SR/du' is parallel (or anti-parallel 

A A A A A A A A A A 

to L and LL:\7M + ML :VL = L-V(L'M) = 0. 

'Wv. *TKM* /V*. /W A*- Ab, 

When ( 36 ) and (39) are added there is a term containing the 

Z A A 

factor [(N/dR/dzi 1 ) + (M-y^/B)] which vanishes because of the relation 
(see reference 7, page 5l) Vp/B - X BR/dz] 1 which holds for these 
(os g) systems. The result for J(z') is 
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^>= J“’<tu'C±)2*Ca*' - z,'B(R)]*L(R)- -T^d u '{[<i> „■ 



u , ,) du " ] _L_ + [<p>u , u-Jdu"] } 


(40) 


C(±)2^(za' ( ~Z 3 ’B)^L' ^-] +0(e a ). 

j(z') can be expressed in terms of the (y 1 , u') variables, 
which have more physical appeal than the z* ones. The reversion 

<wv 

is not too difficult if liberal use is made of the theorem regarding 
f and g type functions. We give only the result: 


J*W'(±)2*{k -MB[R(a, 6, M, K, u")]}*£(R)* ^T 


T ^ B g M? ^ J du”[<ex*(y'> u") - <«> g (y ' , u")j + 0(e 2 ) , 


(41) 


Definitions and equations needed to interpret Eq . (4l) are: 
y': a vector with components a, g, M, and K. 

u 1 : defined in Eq. (20). 

a and g: functions of position, constant on a field line, such that 

A = <*g- 

M: the magnetic moment constant through order e, given by Eq. (l8). 

K: the particle energy, (£ 2 + H 2 )/2. 

T: the period of oscillation, defined below Eq. (20). 

&(y' , o') and j3(y\ u') : the rate of change of a and 3 due to the 

drifts, and given by Eq. (22) and (14a). 
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<a>: the rate of change of <y averaged over a longitudinal 

l 

oscillation and given by > U 1 ). 

o 

<g>: similar to <o>« Thus <a> and <g> are functions of y 1 in 

Eq. (41). 

E, H, J?: the guiding-center variables, defined by Eq, (13). 

jR(a> , M, K, u'): defined by Eq. (26). 

a 

L: the unit vector B/B. 

e: mass to charge ratio m/e. 

The first integral (41) contains e terms in M, which is Mo 4 eMi *, 
for a given M and K the integral is a function of only a and g -- 
i.e., of the field line on which the guiding-center is located, 
and not of u* > which gives the position along the field line. 

The integral is taken between zeros of K — MB(R), even though 
the guiding-center never actually moves along the line (os g) between 
these "virtual 1 ' mirror points, C and D in Fig. 6. The guiding-center 
will eventually be reflected on some field line, but it generally 
will not be the line CD. 

The term of J proportional to T is of order e, since a and g 
are of order e. It depends on u 1 through the upper limit of the 
integral, and vanishes at u 1 r 0, and 1. In fact because eH) 

vanishes when v^ =0, the entire eJi vanishes at the mirror points. 
Thus the guiding-center trajectory will look somewhat as in Fig. 6, 
oscillating about the surface of constant Jo, but always coinciding 
with it at the mirror points. In the figure, the dotted part of the 



Figure 6: The guiding center oscillates about the surface of 

constant Jq. 
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trajectory lies below the Jo surface. 

As shown by Eq. (21) the radical ±2^(K -MB)^ in the first 
term of J/m is not the parallel guiding-center velocity through 
order e. However the integral over a complete period of the 
difference does vanish; the order e term on the right side of (21) 
is proportional v^ and therefore has opposite signs in the ranges 
u' = 0 to J and ^ to 1, whereas iL'dR/du 1 has the same sign. 

Thus the integral of the order e difference vanishes by symmetry. 
For a similar reason, even the term of that is porportional 
to H would contribute nothing to the integral, leaving only the 
sML:VLH 2 Z/B 2 term. 


IV: DIRECT DERIVATION 

Verification that d(J/m)/dt = 0(e 2 ) will serve as a check 
on the result and at the same time suggest a direct derivation of 
J that shortcuts much of the present work. We have 


d J . d ,J N 
dt m " ^da m 


£ 


ii) +°U 2 -> 


(42) 


So we need d(j/m)/da and d(j/m)/d{3 correct only to zero order, 
since a and 0 are of order e. In reference (3) it is shown that 


iL A 

da m 



A Ifi. 

d|3 m 



(43) 
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Furthermore, 0* is 1/T and S(J/m)/^o , equals (T 2 /e) (<j3>& - <a>g) . 
Putting it all together gives what we want: d(J/m)^lt - 0(e 2 ). 

The direct derivation goes as follows: start with the time 

derivative of the lowest order J, 


dt m da m ° dg m 0 3Mo m * 


( 44 ) 


and by (43) convert this to 


^ +Mo + 0(eJ) 


where Mo has been replaced by - cMt . The period T, which depends 
on time has been placed inside the d/d t because its time derivative 
is proportional to the order e drifts, making the error of order e . 
The term containing Mi may be written as 




d 

dM 0 


Jo 

m 



d( Jp/m) -, 
dM 0 J 


- eM, 


_d_ d(Jo/m ) 
dt dM 0 ‘ 


Because Jo is an integral along a field line, so is d(Jo/m) /dMo , 
which is proportional to <B>, the average magnetic field over an 
oscillation. Therefore d(Jo/m)/dMo is changed only by the order e 
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drifts so that eMi “ is °f °(e Z ) 

Transposing all terms of (45) to the 
over time gives 


and can be dropped, 
left side and integrating 


•iDL 


m 


+ eMi 


d(Jo /m ) 

dMo 


d 5' 

P v-'n 


dt[ < p> a — = constant 


or 

j^du' (±)2^[K — (Mo + ~ ^dt[<(3>a — <a>p] - constant 


(46) 


At this point we would have the- result, except for the fact that 
the integral is over time in ( 46 ) and so is not a local quantity. 

This objection can be circumvented by realizing that the correction 
eJi to Jo should be small, which means that the guiding center 
should be distant only order e from the surface on which Jo is 
constant. The time integral from A to F in Fig. 6 can be replaced 
by the integral from C to F along the instantaneous field line, and 
finally we have (41 ) . The u'-dependent term is just the effect of 
the drifts, and the total eJi arises from this and from the correction 
of the lowest-order magnetic moment. 

The direct derivation also reveals that the closed field-line 
case must yield the same expression for J 1 as the oscillatory case, 
which has been considered to this point. In the closed field-line 
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case j the guiding center traverses the line always in the same 
sense 9 parallel or antiparallel to , but with varying speed, 
depending on the magnitude of J3. At the same time it drifts 
slowly at right angles to the line. There is nothing in the 

direct derivation that relies on the oscillatory motion, thus 

we may conclude that the form of the expression ( 41 ) for J will 

be unchanged; u' may be taken as 0 and 1 at an arbitrary point 

on the field line, and only one of the plus -or-minus signs will 
be needed. 

There is however a difference between the oscillatory and 

closed field line cases - namely, that the order e difference 

1 i 

between the integral of v^ and of ± 2 * (K — MB)^ no longer vanishes, 
and the first integral on the right side of ( 4 l) may not be replaced 
by j>v ds, nor may the part of eMi that is proportional to H be 
omitted 

The direct derivation also raises the question of why the time 
derivative of the magnetic moment (as in Eq. 44) did not need to 
be considered in earlier proofs (references 3 and 7 ) of the con- 
servation of Jo. The answer depends on whether Jo is defined with 

(K -MoB)* or with (K —MB)^, which includes some terms of order e - 

* 

If the former, the proofs should indeed consider the effect of Mo. 

Since the proofs show only that <d(Jo/m)/dt> - 0(e 2 ) -- i.e., that 
Jo is only conserved on the average, it is only necessary to show 

• r£ 

that <Mo> vanishes. This is easy to do: <Mo> = — <eM!> - <e/T) \ dtH^ ~ 

(e/T)AM 1 , where AMi is the change in Mt between the time the guid- 
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ing center leaves one mirror and (to lowest order) returns to 
it again. But from (l8) this change is zero because M -j is 
zero when H is zero. If on the other hand Jo is defined with M 
instead of Mo, previous proofs are valid. In practice , defining 
Jo with M is preferable because of the use to which Jo is usually 
put. The invariance of J 0 is often used to determine the surface 
on which the guiding-center moves on the average. In a numerical 
calculation to find the surface one picks a number representing 
the magnetic moment and holds this number constant, so that in 
effect one is using Mo 4- eMi . 

The above discussion leads to the final subject to be treated — 
namely, the order to which the second invariant is conserved. 

Since Jo is an integral along a field line, the parallel velocity 
does not change it, and therefore d(Jo/m)dt = 0(e): Jo is trivially 
conserved to lowest order and is affected by the drifts, which are 
of order e. The usual proof of the invariance of Jo (ref. 3) does 
not show that the e term is zero, but only that its average vanishes 
<d(Jo/m)/dt> = 0(e 2 ). So far it would not matter whether the com- 
plete set of equations (14) were used or the hybrid set (17). The 
drifts are the same for each, and these are what are involved in 
the proofs that Jo is conserved on the average. In the present 
paper we have shown that (l/m) d(Jo 4* eJi)/dt = 0(e 2 ) without 
any averaging involved, but the full set of equations ( 14 ) must be 
used in order that dM/dt = 0(e 2 ). 
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